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During the past few decades, wave energy has received significant attention for harnessing ocean energy. 
Industry has proposed many technologies and, based on their working principle, these technologies 
generally can be categorized into oscillating water columns, point absorbers, overtopping systems, and 
bottom-hinged systems. In particular, many researchers have focused on modeling the point absorber, 
which is thought to be the most cost-efficient technology to extract wave energy. To model such devices, 
several modeling methods have been used such as analytical methods, boundary integral equation meth- 
ods and Navier-Stokes equation methods. The first two are generally combined with the use of empirical 
solution to represent the viscous damping effect, while the last one is directly included in the solution. To 
assist the development of wave energy conversion (WEC) technologies, this paper extensively reviews 
the methods for modeling point absorbers. 
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1. Introduction 


While development of modern wave energy converter dates 
back to 1799 [1], the technology did not receive worldwide atten- 
tion until the 1970s when an oil crisis occurred and Stephen 
Salter published a notable paper about the technology in Nature 
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in 1974 [2]. In the early 1980s, after a significant drop in oil prices, 
technical setbacks and a general lack of confidence, progress slowed 
in the development of wave energy devices as a commercial source 
of electrical power. In the late 1990s, awareness of the depletion 
of traditional energy resources and the environmental impacts of 
the large utilization of fossil fuels significantly increased, thereby 
facilitating the development of green energy resources. The devel- 
opment of the wave energy technology grew rapidly, particularly 
in oceanic countries such as Ireland, Denmark, Portugal, the United 
Kingdom, and the United States. Several pre-commercial devices 
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were deployed. For example, a United States company, Ocean 
Power Technology, deployed one of their 150 kW wave energy con- 
version (WEC) systems in Scotland in 2011 [3]. An Irish company, 
Wavebob, tested a one-quarter scale model in Galway Bay, Ireland, 
in 2006 [4]. In Denmark, the half scale 600 kW Wave Star energy 
system was deployed at Hanstholm in 2009 [5], and a quarter 
and a half size model Wave Dragon was tested at Nissum Bred- 
ning in 2003 [6]. Furthermore, international organizations, such 
as the International Energy Agency and the International Elec- 
trotechnical Commission, are heavily involved in the development 
of wave energy devices. In 2001, the International Energy Agency 


Power conversion 
modules 


Duck rotates 
with nodding 
motion as 
wave passes. 


established an Ocean Energy System Implementation Agreement to 
facilitate the coordination of ocean energy studies between coun- 
tries [7]. In 2007, the International Electrotechnical Commission 
established an Ocean Energy Technical Committee to develop ocean 
energy standards [8]. 


1.1. Device design 
To date, there are more than one hundred prototypes of various 


WEC systems [7]. They can be sorted into a few types. According to 
their energy conversion principles, WEC systems include oscillat- 
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Fig. 1. An illustration of the floating-pitching devices: (top) multiple degree of freedom and (bottom) single degree of freedom. Illustration by Al Hicks, NREL. 
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ing water column, overtopping system, pitch device, membrane, 
and point absorber. Most of these can be both bottom-mounted 
and floating. A brief introduction on the working principle of 
different types of WEC system is given in this section, but the focus 
of this paper is on point absorbers. 

The floating-pitching device, whether it is composed of a sin- 
gle body or a number of connecting bodies, has rotational freedom 
(Fig. 1). The device converts wave energy from its pitching motion. 
The principal axis for floating-pitching devices is either perpen- 
dicular (terminator) or parallel (attenuator) to the wave direction. 
Among various developments since the 1970s, Salter proposed 
one of the most significant inventions in 1974 [2], which became 
the well-known Salter’s Duck (also referred to as the Edinburgh 
Duck). The Salter’s Duck and PS Frog Mk 5 [9] are two examples 
of the terminator, and McCabe Wave Pump [10] and Pelamis [11] 
are the examples of the attenuator. A wave farm consisting of 
750kW Pelamis devices was tested in northern Portugal in 2008 
and another farm was deployed in Orkney, Scotland, in 2010. 

The bottom-hinged pitching system includes a paddle, or a flap, 
that is connected to a hinge deflector on the seabed; the top of the 
device is generally above the free surface (Fig. 2). It is sometimes 
called a wave surge converter, because it converts wave energy 
from the horizontal movement of the water particles. Unlike the 
floating-pitching device, one end of the bottom-hinged device is 
fixed. That is, it shares a similar working principle as the float- 
ing terminators, e.g., Salter’s Duck in Fig. 1(bottom). Examples of 
the bottom-hinged pitching system include the swing mace, also 
proposed by Salter in the 1990s [12], and the Aquamarine Power 
Oyster [13], which was connected to the electrical grid and tested 
in Scotland in 2009 [14]. 

The oscillating water column includes a special chamber with 
a bidirectional turbine inside (Fig. 3). One end of the chamber 


Sea water 
piston 


has an inlet that allows the incident wave to enter and the other 
end contains the turbine. The device converts wave power by 
utilizing the wave elevation to compress or decompress the air 
in a chamber. The compressed air goes through a bi-directional 
turbine. The turbine is a Wells turbine conceived by Professor Alan 
Wells of Queen’s University, in Belfast, in late 1970. A wide variety 
of oscillating water column devices include the LIMPET (shoreline 
system) [15], and Sakata (integrated into breakwater) [16]. 

The overtopping device includes a large structure that embraces 
the incident wave and an outlet with turbines inside the large struc- 
ture (Fig. 4). The device converts wave power by utilizing the wave 
overtopping phenomenon to let the water fall through the outlet 
of the designed structure. When the water falls through the out- 
let, it passes one or more turbines similar to a traditional hydro 
dam; the potential energy is converted into electric power. The 
design involves both kinematic energy and potential energy in the 
conversion process. Examples include the Tapchan (shoreline sys- 
tem) [17], the Wave Dragon (offshore floating system) [18], and 
the Seawave Slot-Cone Generator (SSG), which was integrated into 
a breakwater [19]. 

The membrane device includes a membrane structure and a 
power conversion system that can be a turbine, piezoelectric, or 
other system (Fig. 5). The device converts wave power by utiliz- 
ing the dynamic pressure change during wave propagation (e.g., 
Lilypad [20]). 

A point absorber generally converts wave energy by utilizing 
the heave motion (Fig. 6). Note that there are few pitching point 
absorbers, such as PS Frog Mk5 [9], and they are considered as 
pitching devices in this paper. Typically, a point absorber is either 
a single-body that generates energy by reacting against a fixed 
seabed frame, or it is a multiple-body device that generates energy 
from the relative motion between the two bodies. There are many 


Fig. 2. An illustration of the bottom-hinged system. Illustration by Al Hicks, NREL. 


Y. Li, Y.-H. Yu / Renewable and Sustainable Energy Reviews 16 (2012) 4352-4364 4355 


wells 
turbine 


air flow 
in-take 


column 


Fig. 3. An illustration of the oscillating water column. Illustration by Al Hicks, NREL. 
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Fig. 4. An illustration of the overtopping device. Illustration by Al Hicks, NREL. 
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Fig. 5. An illustration of the membrane device. Illustration by Al Hicks, NREL. 


popular devices such as the OPT PowerBuoy [16], Wavebob [21], 
and Inter Project Service buoy [22-24]. Since the point absorber is 
the focus of the paper, its principles are described in greater detail 
in Section 2. 


1.2. Objective of this paper 


Several literature reviews of WEC devices have been published 
providing information about various aspects of the technology. For 
example, Evans reviewed the analytical derivation of the absorber’s 
motion [25]. Falcdo provide a comprehensive overview of the his- 
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Fig. 6. An illustration of the point absorber. Illustration by Al Hicks, NREL. 


tory and status of the development of WEC systems [26]. Sarmento 
discussed the non-technical aspects of the technology during its 
commercialization path [27]. The U.S. Department of Energy (DOE) 
developed a database to document various technology character- 
istics [28]. Some concise reviews reporting interim progress were 
written by several researchers, such as Falnes [29], on theoretical 
limits [29], and Khan et al. [30], on interconnection issues. There 
also are introductory textbooks on this subject such as Wave Energy 
Conversion by Brooke [31], Ocean Wave Energy by Cruz [32], Ocean 
Wave Energy Conversion by McCormick [33], and Ocean Waves and 
Oscillating Systems by Falnes [34]. 

However, to date, there has been no systematic review of mod- 
eling methods. The rapid development has not led the industry toa 
commercial stage yet. From a pre-commercial stage to a full-scale 
commercial stage, a good understanding of device performance and 
reliability is required. This requirement has not been met yet. In 
addition to these non-technology barriers, it is believed that a lack 
of understanding and quantification of the device behavior is one 
of the main reasons industry development has been delayed. 

During the development of wave energy technology, many 
methods for modeling WECs were developed during the technol- 
ogy’s evolution. To understand point absorber behavior, precise 
modeling tools are needed to simulate the dynamics of the device, 
along with a well-defined combination of modeling methods. Tech- 
nology developers, engineers, and researchers need a guideline 
of which methods to employ for a specific purpose. To facilitate 
the development of the technology and help to model the device 
behavior more cost-effectively, we conducted a systematic review 
of various methods for modeling WECs. The reviews are summa- 
rized in this paper. Since there are a great number of device types, 
we will focus on the point absorber as an example. Additionally, we 
focus on the hydrodynamic modeling methods in this paper. Other 
modeling methods, such as electrical or control modeling, are not 
discussed, although they are important as well. 

Specifically, this paper summarizes the numerical methods used 
for simulating point absorbers. After presenting the working prin- 
ciples of the point absorber, the paper details the existing modeling 
work that uses analytical approaches, empirical methods, and 
numerical methods. At the end, we compare the advantages and 
disadvantages of different hydrodynamic modeling methods. 
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Fig. 7. The translation and rotation of the body (6 degrees of freedom). 


2. Fundamentals of point absorbers 


The simulation of point absorbers is a wave-body interaction 
topic. It requires knowledge of body dynamics and wave theory. 
Before we review the modeling methods for point absorbers, we 
review the fundamentals of body dynamics and wave theory. 


2.1. Point absorber dynamics 


The motion of a point absorber in six degrees of freedom is plot- 
ted in Fig. 7. The dynamics of the device can be expressed by the 
following equations, 


mpat = Frag + Fres + Fpro + Fext, (1) 
Lago + Q x IQ = Mra + Mres + Mpro + Mext, 


where mp is the mass of the body, a; is the acceleration vector 
for the translation, Q and ag are the angular velocity and accel- 
eration vectors, I, is the moment of inertia tensor at the center 
of gravity, F and M are the resulting force and moment acting on 
the body. The subscripts ‘hd’ and ‘res’ represent the hydrodynamic 
excitation component, the restoring component, respectively. The 
subscripts ‘PTO’ and ‘ext’ represent the power take off (PTO) forces 
and the addition constraints, such as mooring line forces. Faltin- 
sen [35] comprehensively reviewed the wave-induced motion and 
loads on offshore structures due to linear, nonlinear, and viscous 
effects. Note that the paper is mainly concerned with the prediction 
of the hydrodynamic components of the system. However, Evans 
[36] and Pizer [37] showed that additional constraints might fur- 
ther reduce the motion response of the device. The effect of the 
mooring constraint depends on its configuration and design. A dis- 
cussion about the influence of the mooring system on the power 
generation performance of a two-body wave energy system was 
presented in Ref. [38]. 

Generally, single-body point absorber converts energy from the 
relative motion between the floating buoy and the seabed. Two- 
body system extracts energy from the relative motion between the 
two oscillating bodies, a float section, and a reaction section. The 
float is only allowed to move freely in one degree of freedom, with 
respect to the reaction section (generally in heave). When the PTO 
component was included, the additional PTO force Fpro was con- 
sidered. A mass-spring—-damper system can be used to represent 
the PTO mechanism, and PTO force becomes 


Fpro = —Cprouret — KZrel, (2) 


where Ze; and Upe; are the relative motion and velocity of the two 
bodies, kis the spring stiffness, Cpro is the power absorption (damp- 
ing) coefficient. The generated power from PTO is defined as 


Pero = Cprou?,,, (3) 


which is proportional to the square of the relative translational 
velocity of the two sections. To optimize the power generation, 
different control strategies can be applied to the system. An intro- 
duction to latching phase-control was given in Ref. [39], and a 
review on the implementation of a spring in phase control, which 
is called reactive phase-control, was described in Ref. [26]. 


2.2. Free surface wave theory 


The free surface waves are a representation of various forces, 
e.g., the wind or a ship, acting on and deforming the fluid surface 
against the action of gravity and surface tension. Naturally, waves 
occur in all sizes and forms. Depending on their magnitude, they 
act on the fluid along with other environmental conditions, such 
as bottom topography, temperature and fluid density. In general, 
waves can be described as linear or nonlinear, regular or irregular, 
unidirectional or omni-directional. The analytical solutions of free 
surface waves were derived based on the potential flow method, 
where the flow is assumed to be incompressible and irrotational. 
The velocity potential @(X,t) therefore, is obtained by solving the 
Laplace equation 


v*¢=0, (4) 


The exact solution of the above equation is very complex. We 
start the review at the linear wave theory, which is obtained using 
simplified assumptions. Linear wave theory, also referred to as 
small amplitude wave theory and airy wave theory [40,41], is the 
simplest solution (first-order approximation) for the flow field. For 
progressive gravity waves of period T, amplitude H, and water depth 
d, linear wave theory assumes H is much smaller than d and wave- 
length A, and the boundary conditions are linearized. The kinematic 
free surface boundary condition (KFSBC) and the dynamic free sur- 
face boundary condition (DFSBC) are then reduced to 


krsec: 9% 2 _o at z=0, 

oz at 

ad (5) 
DFSBC : -gr t8150 at z=0, 


where gis gravity. Both the linearized KFSBC and DFSBC are defined 
at the mean free surface. The flow properties of linear waves can 
be expressed as 


Wave elevation: n= F cos[k(x —ct)], (6) 
; a gH cosh[k(z+d)]_. 
Velocity potential: g= 20 cosh(kd) sin[k(x — ct)], (7) 
> é r 2 w2 E 
Dispersion relation: cf = g7 peanhtkd), (8) 


xH cosh[k(z + d)] 
T sinh(kd) 


Horizontal flow velocity: u= cos[k(x — ct)], 


(9) 


xH sinh[k(z + d)] 
T sinh(kd) 


Vertical flow velocity: w= sin[k(x — ct)], 


(10) 


where x and z are the horizontal and vertical positions, k (=27r/A) 
is the wave number, c (=A/T) is the wave speed, and œ is 
the wave frequency. Depending on the water depth, the waves 
can be categorized as shallow water waves when d/A < 1/20, as 
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Fig. 8. The applicability of various wave theories [45,47]. 


intermediate-depth waves when 1/20 <d/A < 1/2, and as deep water 
waves when d/A > 1/2. 

When the wave amplitude is not small, the waves become non- 
linear. The linearized KFSBC and DFSBC are not satisfied, and the 
nonlinear KFSBC and DFSBC can be written as 


at z= n(x, t), 


P 1|/(ap\* /æÆ\ | a 
prse: 243 B (2) | i 
where C(t) is constant. 

The nonlinear solution can only be obtained by solving Eqs. 
(5)-(8). The nonlinear higher order solution is generally the super- 
position of additional components, which are at higher wave 
frequencies, on the fundamental linear wave-theory terms. There 
have been numerous nonlinear higher order wave theories devel- 
oped since the 1800s. The Stokes theory [42,43] is the most well 
known and studied. An extension of the fifth-order Stokes waves 
was presented by Skjelbreia and Hendrickson [44], and it has been 
widely used in engineering applications. The wave profile and the 
velocity potential of the fifth-order Stokes waves are given as 


(11) 


C(t) at z=n(x,t), 


5 
Wave elevation: kn= ST An cos[nk(x — ct)], (12) 


n=1 


5 
Velocity potential : z = SoG cosh[nk(z + d)] sin[nk(x — ct)].(13) 
n=1 


The applicability of various wave theories was studied in 
Refs. [45,46], and the figure plotted by Le Méhauté is shown in 
Fig. 8. It illustrates the validity limitation of these approximation 
approaches for different wave conditions. A review of the introduc- 
tion of wave theories and the applications to wave load prediction 
on offshore structures is provided by Sarpkaya [47], and a more 
comprehensive review of the theories and the applications of linear 
and nonlinear waves is given in Ref. [48]. 

In the real world, waves are not monochromatic. There are 
a large variety of waves with different frequencies, phases, and 
amplitudes. For an adequate description of the free surface, a large 
number of waves must be superimposed together. Usually, we 
use spectra to describe them. Furthermore, the behavior of waves 


is rather random. In this situation, the waves are called irregular 
waves. The random wave field can be approximated by using an 
infinite sum of sinusoidal components propagating with various 
wavelengths. Pierson—Moskowitz spectrum [49] and Jonswap 
spectrum [50] are the two most widely used spectra profiles. Refer 
to textbooks such as [48,51] for wave theory and [35,47,52] for 
wave-body interactions. 


3. Analytical maximum power extraction studies 


Analytical methods can be sorted, with or without detailed 
descriptions, by WEC systems. The first analytical study of a point 
absorber’s power output focused on maximum efficiency, without 
the detailed description of the device, in 1975 [53,54]. The lin- 
earized equation of motion for the point absorber device in heave 
is 


(m+ my )Z + Rrż + SZ = Fpro + Fext, (14) 


where m is the mass of the body, m, is the hydrodynamic coef- 
ficient of added mass, R, is the radiation damping coefficient, S 
is stiffness of hydrostatic restoring force, Fpro is the vertical force 
due to power-take-off mechanism, and Fext is excitation force. The 
maximum energy absorption Pmax for different incident wave fre- 
quencies can be reached when the oscillating body system is at 
resonance. Based on the linear deep-water wave assumptions, the 
time averaged optimal power absorption at resonance is 


1 
Pmax = ap, Feel (15) 


An absorption width Lmax is defined as 
Pmax _ À 


Lmax = J = Jm’ (16) 
2 2 

_ pg TH 

J= a (17) 


where J is the wave energy flux for linearized deep water waves. 
Note that Evans [55], Mei [56] and Newman [57] also independently 
derived the same result in Eq. (16) using similar approaches in 1976. 
The maximum power that can be absorbed from waves and the 
absorption width is derived based on linear wave theory and con- 
sideration of wave radiation and diffraction. Budal and Falnes [58] 
proposed another upper limit (also referred to as Budal’s upper 
bound) for wave power that is absorbed by a given submerged body 
volume V 
P, | woegH 
Vv ar? 
where the diffraction effect is neglected, and the size of the body 
is assumed to be small compared with the wave length. The two 
curves represent the theoretical prediction of the maximum wave 
energy that can be captured by a semi-submerged heaving body 
with optimum heaving amplitude. The power output for a practical 
device is below the two upper limits, and the maximum power that 
can be absorbed is about half of Pmax, when a phase control method 
is applied [29]. 


(18) 


4. Modeling methods 


The dynamic response of point absorbers can be calculated 
by solving Eq. (1), and power generation can be predicted using 
Eq. (3). Theoretically, the hydrodynamic forces that affect the 
point absorber motion can be obtained analytically! empirically 


1 The analytical method can be derived by an equivalent linearization of the 
empirical calculation. Therefore, the former can be viewed as a special case of the 
latter. 
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or numerically, depending on the modeling purpose, the opera- 
tional condition, and the detailed design. These modeling methods 
evolved from modeling ship hydrodynamics and offshore floating 
structures. Therefore, we apply the methods for naval architecture 
and offshore engineering when we review modeling methods for 
predicting the hydrodynamics of point absorber systems. 


4.1. Analytical methods 


Acommon approach to determine the hydrodynamic force ana- 
lytically is to presume that the force is the sum of wave diffraction 
and wave radiation components, 


Fha = Fr + Fp, (19) 


where the subscripts ‘R’ and ‘D’ represent the radiation and diffrac- 
tion terms, respectively. The radiation components include an 
added-mass component and a radiation damping term. 


where A and B are the hydrodynamic added-mass and damping 
coefficients. The wave excitation force, induced by wave diffrac- 
tion, includes a diffraction component Fpr, due the present of the 
floating body, and a Froude-Krylov force component Frx, generated 
by undisturbed incident waves. 


Fp = Ode Kr + OD) ads = Fey + Fpr, (21) 
sw t oe o 


where Sw is the immersed surface on the body, n is the unit nor- 
mal vector, ġ; is the incident wave velocity potential and @p is the 
diffraction velocity potential. 

The velocity potential are then solved analytically, and most 
of the analytical solutions were derived through the use of linear 
wave theory. Many analytical studies can be found in the literature. 
For examples, Newman [59] calculated the excitation force and 
moment though the use of Haskind relation. Yeung [60] proposed 
an analytical method for calculating the added mass and damp- 
ing coefficients for offshore engineering applications. The paper 
suggests a method to calculate the coefficient using the Eigen func- 
tions. In addition, an analytical derivation of the hydrodynamic 
coefficients of a floating hemisphere undergoing forced periodic 
oscillations was presented in Ref. [61]. The analytical methods have 
also been widely used in ocean engineering problems, especially for 
an array of cylinders [62,63]. 

In reality, the viscous damping effect is important for modeling 
the hydrodynamics of point absorbers, and neglecting the viscous 
drag will result in an unrealistic large prediction on the response of 
the device and on power extraction performance, particularly near 
the resonance. Discussion on energy loss due to viscous damping 
can be found in [34,64]. A common approach to include the viscous 
drag effect in the modeling is to add a quadratic damping term to the 
non-viscous hydrodynamic forces. More details on the calculation 
of viscous damping force are discussed in Section 4.3. 

In the past decade, many researchers have used the analyt- 
ical method to calculate the power output of point absorbers 
[26,65-67]. In addition to a single unit, this method has been used 
to study arrays of WEC systems. For example, through the use of 
linear theory and multiple-scale approximation, Garnaud and Mei 
[68] presented an analytical solution to the wave energy extrac- 
tion using infinite strips of buoys and a circular array. The study 
showed that an array of heaving wave-absorbing buoys could have 
a broader range of wave frequencies than a single large buoy, with 
good energy-absorbing efficiency. In general, analytical methods 
are very efficient for providing a quick performance estimation of 
point absorbers with simple geometries. 


4.2. Boundary integral equation methods 


Boundary integral equation method (BIEM), also referred to as 
boundary element method, is an advanced potential flow method 
that can be used for handling more complicated geometries. The 
governing equation, Eq. (4), is formulated in a boundary integral 
equation form, which is derived by introducing a Green’s function 
G(p.q), where p and q represent the field point and the source point, 
respectively. The boundary integral equation obtained through the 
use of Green’s third identity gives 


a(p)¢(p) + i P(P)Gn(Pp, q) — On(P)G(p, q) ds = 0, (22) 
S 


where S indicates the boundaries of the fluid domain, a(p) is the 
internal angle formed at the boundaries, Gn = vG-n and @n=V¢@-n, 
with n being the normal vector at q. The boundary surfaces are 
discretized into panels. By using Dirichlet- and Neumann-type 
boundary conditions, the potential flow field can be obtained by 
solving the resulting system of linear equations numerically. 

After the potential flow field is obtained, the pressure on the 
body surface is evaluated using Bernoulli’s equation, 


p=-p($e+ 5V0-Vo+ez). (23) 


The forces F and moment M acting on the body are calculated 
by integrating the pressure on the immersed body surface, 


p= pn ds; M= pn x rds, (24) 
Sg Sg 


where r is the distance vector from the center of rotation. The 
viscous damping coefficient can be estimated based on the empir- 
ical solution, following the same approach used in the analytical 
method. 

Regarding the simulation of the interaction between waves and 
floating bodies, two types of BIEM approaches have been used 
for modeling. One follows a weakly nonlinear theory that uses a 
perturbation expansion for the solution (the theory of wave radi- 
ation and diffraction) and specifies the boundary conditions at the 
mean free surface and body-surface through the use of Taylor series 
expansion. The other method uses a fully nonlinear time-domain 
approach. 

The weakly nonlinear approach has the advantage of having the 
same coefficient matrix for the system of equations solved at every 
time step. One can solve the problem in the time domain [69] or 
in the frequency domain, which is the method used in the panel 
code WAMIT [70]. In the frequency domain calculation, radiation, 
diffraction, and excitation forces are calculated through the use 
of the linear superposition principle, which is limited to smaller 
amplitude waves (weak nonlinear waves) [71]. The complexity of 
the method increases significantly for higher order schemes. Exam- 
ples of studies that used this weak nonlinear approach for wave 
radiation and diffraction problems can be found in Refs. [72-78]. 
A discussion on the use of linear and second-order approaches for 
predicting the wave excitation of floating bodies was given in Ref. 
[79]. 

The fully nonlinear time domain approach for wave analy- 
sis was first proposed by Longuet-Higgins and Cokelet, using a 
Mixed Eulerian-Lagrangian (MEL)-type method for modeling sur- 
face waves [80,81]. Fully nonlinear boundary conditions were 
applied on the instantaneous water free surface and body surface, 
and the free surface was computed using a tracking method. With 
BIEM, the computational domain is discretized along the domain 
boundaries, including the free surface and the body surface. There- 
fore, the coefficient matrix for the system of equations should 
be updated at every time step, and the method is only solved in 
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time domain. A high order Runge-Kutta scheme is often used to 
update the time varying free surface and body surface. The non- 
linear time-domain approach is more accurate for modeling the 
highly nonlinear interaction between waves and floating bodies. 
Using this approach, researchers successfully predicted the nonlin- 
ear irregular wave radiation [82] and diffraction [83] of vertical 
cylinders, as well as the ventilating entry of a surface-piercing 
hydrofoil [84]. 

Additionally, modeling floating bodies in numerical wave tank 
requires a boundary that generates waves and absorbs reflected 
waves at the same time, and a boundary that absorbs outgo- 
ing waves. Over the years, many efforts have been made to 
develop a numerical method to simulate these types of bound- 
aries, without creating additional numerical disturbances. A review 
on unbounded domains, with various types of artificial boundary 
conditions and absorbing methods, was given in [85]. Most of the 
studies between the 1970s and the 1990s were based on potential 
flow theory. In particular, a sponge-layer method was proposed by 
Israeli and Orszag [86], in which an artificial damping layer (zone) 
was implemented to attenuate the outgoing wave before it reached 
the outer boundary. Similar approaches also were proposed in Refs. 
[87,88]. 

The frequency domain BIEM is the most commonly used BIEM 
method for modeling the hydrodynamics and the power genera- 
tion performance of point absorber systems. The method is used 
for calculating the non-viscous hydrodynamic forces described in 
the analytical method. McCabe et al. [89] demonstrated an exam- 
ple of modeling the responses of axisymmetric bodies in irregular 
waves using WAMIT. The method was combined with a Laplace 
transfer-function formulation to model the dynamics of the floating 
body in time domain. Given that the BIEM can be used for modeling 
more complex geometry, the method has been widely used in the 
studies of not only point absorber systems but also a wide variety 
of WEC systems. For examples, the method was used to study the 
power generation performance of floating pitching devices [9,90], a 
modified version of the Salter’s Duck by Cruz and Salter [91], anda 
bottom-hinged surge-pitching plate [92]. Moreover, a recent study 
utilized the frequency domain BIEM to perform a series of studies 
on the power extraction performance of eight selected WEC sys- 
tems, including five heaving absorber systems with different design 
principles [64]. The study also analyzed the annual power gener- 
ation performance and the capture efficiency at five different site 
locations. Moreover, several researchers have used the method in 
the study of control and tuning strategies for heaving absorbers 
[93-96]. For the analysis of array systems for point absorbers, the 
interaction effect could be considered through the use of a multiple- 
scattering method [67,97,98]. In particular, Vicente et al. [98] used 
the method to analyze the dynamics of the array system with inter- 
body connections and bottom slack-mooring lines. 


4.3. Viscous drag calculation for potential flow solution 


Generally, the additional viscous damping term follows the drag 
term in the Morison’s equation. The Morison’s equation was pro- 
posed by four authors, J.R. Morison, M.P. O’Brien, J.W. Johnson and 
S.A. Schaaf, at the University of California, in 1950 [99], and the drag 
term was given as 


1 
Fa = — 5 PCaAc(u — Up)|U— up| (25) 


where C, is the quadratic drag (viscous damping) coefficient, p indi- 
cates the density of sea water, Ac is the characteristic areas, u is 
the flow velocity, and u, is the velocity of the body. Generally, the 
value of viscous damping coefficient depends on the body geom- 
etry, the Reynolds number and the Keulegan—Carpenter number. 


More details of the use of quadratic drag approach and the appro- 
priate values for the empirical coefficient can be found in [47,100]. 

In addition, the empirical viscous damping coefficient data are 
limited to existing simple geometry. For practical point absorber 
geometry, the hydrodynamic forces may have to be evaluated by 
conducting additional wave tank tests [101] or prescribed motion 
CFD simulations [102,103] if existing data is not available. 


4.4. Navier-Stokes equation methods 


The viscous effects of boundary layer separation, turbulence, 
wave breaking, and overtopping are important to the prediction of 
hydrodynamic forces of the devices. The potential flow method can- 
not capture them. Therefore, modelers often turn to a fully viscous 
solution by implementing the Navier-Stokes equation methods 
(NSEM). The velocity vector U and pressure p of the incompress- 
ible flow field are obtained by solving the continuity equation and 
the momentum equations, which are given as 


V-U=0, 


p(U; + U- VU) = (26) 


Vp+F,+V-T, 


where U; is the time derivative of U, F, is the body force vector (such 
as gravity), and T indicates the stress tensor. The computational 
domain can be discretized using a finite difference/volume/element 
method, and the resulting system of linear equations can be solved 
directly, or preferably, through iterations using an appropriate lin- 
ear equation system solver. A well known Semi-Implicit Method for 
Pressure-Linked Equations (SIMPLE) algorithm [104] and its vari- 
ants have been widely used for solving Navier-Stokes equations. 
More details on the numerical discretization methods, as well as 
the physical phenomena, are given in Refs. [105,106]. The force 
and moment calculation is similar to that applied in the potential 
flow theory, except that the frictional force also is included in the 
calculation. 

Basically, using the Navier-Stokes equation-based method for 
modeling the hydrodynamics of floating bodies in waves involves 
the calculation of the free surface in a numerical wave tank and 
the simulation of the turbulent flow. In general, two types of 
approaches have been applied for predicting the free surface. One 
is the surface tracking method, which treats the free surface as 
a sharp boundary and updates it with time [107]. However, the 
method does not have the ability to model wave breaking and 
overtopping. The other is the interface-capturing method, where 
the simulation is performed on a grid that includes both air and 
water phases. Marker-and-Cell (MAC) method [108], Volume Of 
Fluid (VOF) method [109], and the level set approach [110,111] 
are the most often used interface-capturing methods. In particu- 
lar, the VOF method is the most widely used for commercial and 
open source computational fluid dynamics (CFD) codes. Review of 
these methods and their application to wave hydrodynamics can 
be found in [106,112]. The artificial damping layer method also is 
widely used in the Navier-Stokes equation-based approaches for 
absorbing the outgoing and reflecting waves in the wave-body 
interaction simulation. For example, Lara et al. [113] combined the 
sponge-layer method, with an internal-wave maker method [114], 
to model the wave generating-absorbing boundary and simulated 
the interaction between waves and permeable structures. 

In general, four types of numerical methods are applied 
for turbulent flow modeling, including direct numerical sim- 
ulations (DNS) [115], large eddy simulations (LES) [116,117], 
Reynolds-averaged Navier-Stokes (RANS) methods, and detached 
eddy simulations (DES). The corresponding computational cost for 
DNS and LES is very high DNS, and the two models are often 
viewed as a research tool in the study of turbulence, particularly 
for most floating body dynamics problems. Therefore, the RANS 


Y. Li, Y.-H. Yu / Renewable and Sustainable Energy Reviews 16 (2012) 4352-4364 4361 


method is the most widely used scheme for modeling turbulent 
flow. The unsteadiness of the turbulence is time averaged and 
the computation of turbulent flow relies on simplified engineering 
approximations. DES is a hybrid method that combines the RANS 
and LES methods. It handles the near wall region, with a RANS- 
type model, and treats the rest of the flow field using a LES-type 
method. Different types of RANS models are designed for different 
turbulent flows. The physical phenomena and the basic assump- 
tions for each turbulence model were comprehensively reviewed 
by Wilcox [118]. 

Using NSEM for modeling floating body dynamics often involves 
the handling of moving boundary and mesh deformation. Sev- 
eral methods have been developed over the years, including: 
(1) boundary-fitted mesh morphing method, which is the most 
common approach and typically involves the use of Arbitrary 
Lagrangian-Eulerian methods [119]; (2) overset method (also 
referred to as Chimera approach), which use a multi-mesh system 
that solve the problem using separately generated but overlapping 
meshes through iterations [120]; (3) immersed boundary method, 
which introduces a body force term in the governing equation to 
represent the boundary surface [121]; (4) Smooth Particle Hydro- 
dynamics (SPH) method [122], which is a mesh free Lagrangian 
method and uses a smoothing kernel to approximate a delta func- 
tion, with many sampling points in the fluid domain. Generally, 
boundary-fitted mesh morphing method and the overset method 
are the most widely used models for solving wave and floating body 
dynamics problems. 

NSEM has often been used for studying the complex nonlinear 
wave and floating body interaction in ocean engineering problems. 
For example, Moctar et al. [123] used a RANS model to analyze the 
wave load of a jack-up platform under freak waves, where the VOF 
method was adopted for the free-surface calculation and was cou- 
pled with a finite element model for the structure’s stress analysis. 
The study demonstrated the significance of using a Navier-Stokes 
type approach, where the effects of wave run-up on the platform 
leg and the impact-related wave loads on the hull were all con- 
sidered. Their study also showed that the difference between the 
shear and overturning moment for the platform under freak waves 
calculated from the Morison’s approach, and those obtained from 
the RANS calculation, can increase up to 25%, especially for larger 
and breaking waves. 

A few studies on the fluid-structure interaction of WEC sys- 
tems were performed using the Navier-Stokes-type approaches. 
A RANS code (COMET) was applied by Agamloh et al. [124] to 
model the dynamics of a single buoy and a double buoy system 
in waves. The PTO mechanism was considered as the damping 
system based on the energy that was captured. Studies on the 
Wave-driven, Resonant, Arcuate-action, Surging Point-Absorber 
(WRASPA) were presented by Bhinder et al. [125,126]. The authors 
performed a series of studies using a RANS model to investigate 


the hydrodynamics of a floating-point absorber [127]. The power 
take off was represented using a mass-spring—damper system, and 
the maximum power that this particular device generated, with a 
constant value spring—damper system, was examined. 


5. Discussion 


In Section 4, we have summarized the existing numerical meth- 
ods for modeling the dynamic response of point absorber in waves 
and for predicting the power generation of the system. Here, we 
provide further discussion and insights on these methods. 

Table 1 illustrates the features of each method for modeling 
point absorber systems. Depending on the purpose of the study, 
each method has its advantages and limitations. The analytical solu- 
tion only can be used for simple geometries, while BIEM and NSEM 
can be used for more complicated and realistic ones. The analytical 
method and the frequency domain BIEM model the system dynam- 
ics in a sequential approach and use the treatments of radiation 
and diffraction problems for estimating the hydrodynamic forces. 
Therefore, the utilization of these two methods is generally limited 
to relatively smaller waves. On the contrary, the time domain BIEM 
and NSEM fully couples the flow field simulation with the solution 
for the equation of motion at each time step. The time domain BIEM 
is feasible for capturing the fully nonlinear interaction between 
waves and the floating body, except for wave breaking, overtopping 
and viscous flow separation. A floating body dynamics CFD (NSEM) 
simulation is needed to fully capture all the nonlinear interaction 
between waves and the device, including viscous flow separation, 
turbulence, wave breaking and overtopping. 

An example for modeling the power generation performance of 
a two-body floating-point absorber system is presented to demon- 
strate the utilization of these methods. The two most commonly 
used methods, a frequency-domain BIEM and a NSEM (CFD), were 
selected. In the frequency domain BIEM, the hydrodynamic forces 
were calculated using WAMIT, and the viscous damping coeffi- 
cient was estimated based on the empirical solution of [47,100]. 
The details of the device geometry and the numerical settings in 
the CFD simulation were described in Ref. [127]. To compare the 
results obtained from the two methods, the computational time for 
each model is shown in Table 2, and the power generation perfor- 
mance is plotted against incident wave periods in Fig. 9. It is found 
that the frequency-domain BIEM is much more cost-effective. How- 
ever, the selection of viscous damping coefficients has a significant 
influence on the frequency-domain BIEM results, particularly near 
resonance. Note that the maximum power generally occurs when 
point absorber systems are close to resonance. The calculation of 
viscous damping is important to predict the power and may lead 
to significant uncertainty in the prediction results. For example, 
Babarit et al. [64] showed that if the drag coefficients were varied 
by a factor of 4, the effect on power generation could be up to 30%. 


Table 1 
Features of hydrodynamic modeling methods. 
Methods Features 
Viscous effects Dynamic coupling Wave/body interaction Accuracy 
complexity 
Hydrodynamics modeling System dynamics modeling 
Analytical Empirical Radiation and Linear Affected by the prediction Limited to waves with 
diffraction theory of viscous damping relatively smaller 
Frequency domain Empirical Radiation and Non-breaking waves Affected by the prediction amplitude 
BIEM diffraction theory of viscous damping 
Time domain BIEM Empirical Fully coupled time Non-breaking waves Affected by the prediction Very good 
domain simulation of viscous damping Fully nonlinear 
NSEM: floating-body Solves Fully coupled time Fully nonlinear and Very good 


Navier-Stokes domain simulation 


equation 


dynamics CFD 


breaking waves 
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Table 2 
Run-time benchmark for modeling point absorber WEC systems. 


Hydrodynamics modeling Discretization 


Number of panels/meshes Wall clock time 


Frequency domain BIEM (WAMIT) & quadratic drag 
NSEM: floating body-dynamics CFD Computational domain 


discretization 


Boundary discretization 


5096 (low-order panels) 
~1 million (finite volume mesh) 


2804s? (30 wave frequencies) 
12h on 64 cores? (one wave frequency) 


a On a 3.33 GHz Intel i5 processor with 8 GB of memory. 


b Each compute node consists of dual socket/quad-core 2.93 GHz Intel Nehalem processor, with 12 GB of memory shared by all 8 cores. 


30 


Power/H? (kW/m?) 


10 12 
Wave period (sec) 


CFD from Yu & Li (2011) 
WAMIT+vis, (0.5xCd) 
WAMIT+vis, (1xCd) 
WAMIT+vis, (2xCd) 


Fig. 9. Power generation prediction from NSEM and BIEM (Cpro = 507 kN s m~'). 


For operational wave conditions, where waves are relatively 
small and linear theory are generally applicable, and the non- 
viscous hydrodynamic forces are dominated by radiation and 
diffraction components. The frequency domain BIEM is the most 
commonly used method for modeling the power generation per- 
formance of point absorber systems. However, the viscous damping 
coefficients need to be carefully selected. The drag coefficient 
depends on the device geometry, the scale and the relative 
velocity between the body and the flow around it. The drag coef- 
ficient becomes much larger when the Reynolds number and 
the Keulegan-Carpenter number are smaller [100]. In addition, 
after the frequency domain hydrodynamic forces are obtained, the 
response of point absorbers can be calculated using a frequency 
domain approach or a time domain approach. The frequency 
domain solution represents the steady-state response due to the 
sinusoidal excitation of regular waves. The time domain approach 
solves the Cummins equations [128] and are often used for mod- 
eling the dynamics of the system in irregular waves, particularly 
for the study of optimal control strategies. The methods and the 
relationship between the parameters in the frequency domain 
approach and those in the time domain approach were well sum- 
marized by Ogilvie [129]. 

For survivable conditions, a time-domain fully nonlinear BIEM 
model or NSEM is more applicable. In reality, point absorbers are 
designed to sustain the extreme wave load during storms. The 
device may subject to waves with a typical 100-year significant 
wave height in the range between 8 m and 13 m. Under these waves, 
breaking and overtopping waves often occur and the effect of the 
nonlinear interaction between waves and the floating body on the 
dynamic response of the system is significant. Only NSEM models 
have the capability of modeling all the complex nonlinear hydrody- 
namics. However, longer computational time is required. The actual 
computational time generally depends on many factors, including 


processor speed, convergence criterion, mesh quality, parallel effi- 
ciency, and the complexity of the interaction between waves and 
the floating body. Thus, one may have to compromise accuracy 
against computational cost. 


6. Conclusion 


Four numerical methods for modeling the point absorber sys- 
tems have been reviewed in this paper; each method has its own 
strength and limitation. For idea analysis, analytical method is suf- 
ficient and effective. For modeling a practical point absorber WEC 
system, frequency domain BIEM and NSEM are the most commonly 
used method depending on the wave condition. In operational 
condition, the frequency domain BIEM, with a carefully selected 
viscous damping, is generally used for power generation perfor- 
mance analysis, design and optimization. They tends to be easier 
to combine with the modeling techniques for other aspects of the 
point absorbers such as the power system. On the other hand, NSEM 
is typically used for detailed analysis or for survivability analy- 
sis, where the nonlinear interaction between waves and the WEC 
system is significant. Time domain fully nonlinear BIEM is more 
appropriate in the transition between these those conditions and 
for developing active control strategies. 

In additional to the point absorbers, there are a wide variety of 
WEC systems such as bottom hinged device and oscillating water 
column. Similar modeling philosophy can be applied to them. How- 
ever, the dominant physical phenomenon for each system varies. 
For examples, additional nonlinear effects on the hydrodynamics 
of the system due to wave breaking and overtopping, can be signif- 
icant for a bottom-hinged pitching plate. The additional constraint 
due to mooring connections can be essential for a floating device 
oscillating water column. The prediction of these additional forces 
and the level of complexity of the interaction between forces and 
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system dynamics will affect the accuracy of the modeling method, 
particularly the frequency domain BIEM, and the overall computa- 
tional time. We will investigate these features in future. 
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